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Abstract. We study a class of discrete-time multi-agent systems modelling opinion dynamics with 
decaying confidence. We consider a network of agents where each agent has an opinion. At each time 
step, the agents exchange their opinion with their neighbors and update it by taking into account 
only the opinions that differ from their own less than some confidence bound. This confidence bound 
is decaying: an agent gives repetitively confidence only to its neighbors that approach sufficiently fast 
its opinion. Essentially, the agents try to reach an agreement with the constraint that it has to be 
approached no slower than a prescribed convergence rate. Under that constraint, global consensus 
may not be achieved and only local agreements may be reached. The agents reaching a local agreement 
form communities inside the network. In this paper, we analyze this opinion dynamics model: we 
show that communities correspond to asymptotically connected component of the network and give an 
algebraic characterization of communities in terms of eigenvalues of the matrix defining the collective 
dynamics. Finally, we apply our opinion dynamics model to address the problem of community 
detection in graphs. We propose a new formulation of the community detection problem based on 
eigenvalues of normalized Laplacian matrix of graphs and show that this problem can be solved using 
our opinion dynamics model. We consider three examples of networks, and compare the communities 
we detect with those obtained by existing algorithms based on modularity optimization. We show 
that our opinion dynamics model not only provides an appealing approach to community detection 
but that it is also effective. 



1. Introduction 

The analysis of multi- agent systems received an increasing interest in the past decades. In such 
systems, a set of agents interact according to simple local rules in order to achieve some global 
coordinated behavior. The most widely studied problem is certainly the consensus or agreement 
problem where each agent in the network maintains a value and repetitively averages its value with 
those of its neighbors, resulting in all the agents in the network reaching asymptotically a common 
value. It is to be noted that the graph of interaction describing the network of agents is generally 
not fixed and may vary in time. Conditions ensuring consensus have been established by various 
authors including |JLM03l IBHOT05L IMor05l iRBOSj (see [OSFMOTj for a survey). More recently, 
there have been several works providing estimations of the rate of convergence towards the consensus 
value |OT09[ I^BOSl IZW09j . 

In this paper, we adopt a different point of view. We consider a discrete-time multi-agent system 
where the agents try to reach an agreement with the constraint that the consensus value must be 
approached no slower than a prescribed convergence rate. Under that constraint, global consensus 
may not be achieved and the agents may only reach local agreement. We call communities the 
subsets of agents reaching a consensus. Our model can be interpreted in terms of opinion dynamics. 
Each agent has an opinion. At each time step, the agent receives the opinions of its neighbors and 
then updates its opinion by taking a weighted average of its opinion and the opinions of its neighbors 
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that are within some confidence range of its own. The confidence ranges are getting smaller at each 
time step: an agent gives repetitively confidence only to the neighbors that approach sufficiently 
fast its own opinion. This can be seen as a model for a negotiation process where an agent expects 
that its neighbors move significantly towards its opinion at each negotiation round in order to keep 
negotiating. Our model can be seen as an extension of the opinion dynamics with bounded confidence 
proposed by Krause in |Kra97| and studied in |HK021 IBHT09] . 

We analyze our opinion dynamics model by first studying the relation between asymptotic agreement 
of a subset of agents and the fact that they are asymptotically connected. We show that under 
suitable assumptions, these are actually equivalent (i.e. communities correspond to asymptotically 
connected component of the network) except for a set of initial opinions of measure 0. We then 
give an algebraic characterization of communities in terms of eigenvalues of the matrix defining the 
collective dynamics. 

Finally, we apply our opinion dynamics model to address the problem of community detection in 
graphs. In the usual sense, communities in a graph are groups of vertices such that the concentration 
of edges inside communities is high with respect to the concentration of edges between communities. 
Given the increasing need of analysis tools for understanding complex networks in social sciences, 
biology, engineering or economics, the community detection problem has attracted a lot of attention 
in the recent years (see the extensive survey [ForlOj ) . The problem of community detection is however 
not rigorously defined mathematically. Some formalizations of this problem have been proposed in 
terms of optimization of quality functions such as modularity |NG04j or partition stability |LDB09| . 
We propose a new formulation of this problem based on eigenvalues of normalized Laplacian matrix 
of graphs and show that this problem can be solved using our opinion dynamics model. We consider 
three examples of networks, and compare the communities that we detect with those obtained by 
the modularity optimization algorithms presented in [New061 IBGLL08] . We show that our opinion 
dynamics model not only provides an appealing approach to community detection but that it is also 
effective. 



2.1. Model Description. We study a discrete-time multi-agent model. We consider a set of n 
agents, V = {1, . . . ,n}. A relation E <^ V x V models the interactions between the agents. We 
assume that the relation is symmetric {{i,j) £ E iS {j,i) G E) and anti-reflexive ((i,i) ^ E). V 
is the set of vertices and E is the set of edges of an undirected graph G = {V,E), describing the 
network of agents. Each agent i gV has an opinion modelled by a real number Xi{t) G M. Initially, 
agent i has an opinion Xj(0) = independent from the opinions of the other agents. Then, at every 
time step, the agents update their opinion by taking a weighted average of its opinion and opinions 
of other agents: 
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with R > and p G (0, 1) model parameters. 

Remark 1. It is noteworthy that the confidence neighborhoods Ni{t) and the coefficients Pij{t) ac- 
tually depend also on the opinions xi{t), . . . ,Xn{t). For the sake of simplicity and in order to reduce 
the length of the equations we keep the notations pij{t) and Ni{t) pointing out just the variation in 
time of these quantities. 



We make the following additional assumptions: 

Assumption 1 (Stochasticity). Fort£N, the coefficients pij{t) satisfy 

(a) pij{t) £ [0, 1], for all i,j £ V. 

(b) E]=iPijit) = 'i-, for allieV. 



This model can be interpreted in terms of opinion dynamics. At each time step t, agent i G V 
receives the opinions of its neighbors in the graph G. If the opinion of i differs from the opinion of 
its neighbor j more than a certain threshold Rp^, then i does not give confidence to j and does not 
take into account the opinion of j when updating its own opinion. The parameter p characterizes the 
confidence decay of the agents. Agent i gives repetitively confidence only to neighbors whose opinion 
converges sufficiently fast to its own opinion. This model can be interpreted in terms of negotiations 
where agent i requires that, at each negotiation round, the opinion of agent j moves significantly 
towards its opinion in order to keep negotiating with j. 

This model is somehow related to the one discussed in |CS771 ICHN86j where agents harden their 
position by increasing over time the weight assigned to their own opinion. In our model, the agents 
implicitely increase also the weights assigned to their neighbors whose opinion converges sufficiently 
fast to their own opinion, by disregarding the opinions of the other agents. As noticed in |CHN86j . 
hardening the agents positions may hamper the agents to reach an asymptotic consensus. This will 
be observed in our model as well. However, the aim in this paper is not to exogenously increase the 
self-confidence of the agents, but to meet a prescribed convergence speed towards the final opinion 
profile. 

Remark 2. We assume in this paper that p G (0, 1). However, let us remark that for p = 1 (there 
is no confidence decay), with a complete graph G (every agent talks with all the other agents), and 
with coefficients Pij{t) given for all j G {i} U Ni{t) by 

Pij{t) = -^,\(^^ w^^h di{t) = 1 

our model would coincide with Krause model of opinion dynamics with bounded confidence |Kra97| 
IHK021 IBHT09] . 



Our first result states that the opinion of each agent converges to some limit value: 

Proposition 1. Under Assumption^ (Stochasticity), for all i £ V, the sequence (xi(t))jgN is con- 
vergent. We denote x* its limit. Furthermore, we have for all t G N, 

(2.4) |^^(t)_^*|<_^pt. 
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Proof. Let z G t G N, we have from (2.1), Assumption [T] and (2.2) 



\x^{t + 1) - Xi{t)\ 



^Pij{t)Xj{t) - Xi{t) 



< 



^Pij{t){xj{t) - x^{t)) 

i=i 



^ Pij{t){Xj{t) - Xi{t)) 
j(LN,{t) 

^ Pij{t)\xj{t) - Xi{t)\ 



Then, it follows from equation (2.3) that 



\Xi{t+l)-Xi{t)\ < P^jit)Rp' 

Finally, Assumption [T] gives for all t G N 

\xi{t + I) - Xi{t)\ < {l-pu{t))Rp' < RpK 

Let t G N, T G N, then 

r-1 r-1 



i{t + t) - Xi{t)\ <Y\xi{t + k + l)-Xi{t + k)\ < J2^P 



t+k 



k=0 



k=0 



Therefore, 
(2.5) 



R 

\xi{t + r) - Xi{t)\ < p\l - p^) < 



R 



-P 



1-p^ - ' ' - l-p' 

which shows, since p G (0, 1), that the sequence {xi{t))t^^ is a Cauchy sequence in M. Therefore, it 
is convergent. Equation (2.4) is obtained from (2.5) by letting r go to +c«. □ 



Remark 3. The convergence of each opinion sequence {xi{t))tiz.^ could have been proved using a result 
from |Lor05| . even for p = 1, with the additional assumption that the non-zero coefficients pij{t) are 
uniformly bounded below by some strictly positive real number. However, the result in |Lor05j does 
not provide an estimation of the convergence rate which is essential in our subsequent discussions. 

The previous proposition allows us to complete the interpretation of our opinion dynamics model. 
The agents try to reach an agreement with the constraint that the consensus value must be ap- 
proached no slower than 0(p*). Under that constraint, global agreement may not be attainable and 
the agents may only reach local agreements. We refer to the sets of agents that asymptotically agree 
as communities. 

Definition 1. Let i,j G V, we say that agents i and j asymptotically agree, denoted i ~* j, if and 
only if X* = x*. 

It is straightforward to verify that ~* is an equivalence relation over V. 
Definition 2. A community C '^V is an element of the quotient set ^ = V/ ^* . 
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Let us remark that the community structure is dependent on the initial distribution of opinions. In 
the following, we shall provide some insight on the structure of these communities. But first, we 
need to introduce some additional notations. 

2.2. Notations and Preliminaries. We define the set of interactions at time t, E{t) C y x y as 

E{t) = G E\ \xi{t) - Xj{t)\ < V} • 

Let us remark that G E[t) if and only if j G Ni[t). The interaction graph at time t is then 

G{t) = {y,E{t)). Let us remark that Remark [l] applies also to E{t) and G{t). 

For a set of agents / C 1/, the subset of edges of G connecting the agents m. I \s Ei = E r\ {I x I). 
Let E' C i?/ be a symmetric relation over /, then the graph G' = (I, E') is called a subgraph of 
G. If / = V , then the graph G' = {V, E') is called a spanning subgraph of G. The set of spanning 
subgraphs of G is denoted S{G). For all t G N, G{t) G S{G). Let us remark that the set S{G) is 
finite: it has 2l^l/^ elements (because we only consider symmetric relations) where \E\ denotes the 
number of elements in E. Given a partition of the agents V = {Ii, . . . , Ip}, we define the set of edges 
Ep = IJj^-p Ej and the spanning subgraph of G, G-p = {V,Ep). Essentially, G-p is the spanning 
subgraph of G obtained by removing all the edges between agents belonging to different elements of 
the partition V. An interesting such graph is the graph of communities G^^ = (y, K^) where: 

E^^ = {{i,j)eE\i^* j}. 

Let G' = {V, E') G 5(G), a path in G' is a finite sequence of edges (11,12), (^2, ^3), • • ■ , (ip, ip+i) such 
that (ik,ik+i) G E' for all /c G {1, . . . ,p}. Two vertices i, j G V are connected in G' if there exists a 
path in G' joining i and j (i.e. ii = i and jp = j). A subset of agents I CV is a connected component 
of G' if for all i, j £ I with i ^ j, i and j are connected in G' and for all i £ I, for all j €z V \ I, 
i and j are not connected in G' . The set of connected components of G' is denoted IC(G'). Let us 
remark that /C(G') is a partition of V. 

We define the vectors of opinions x(t) = (xi(t), . . . , x„(t))^ and of initial opinions x^ = (x^, . . . , x^)~^ . 
The dynamics of the vector of opinions is then given by 

x(t + l) = P(t)x(t) 

where P(t) is the row stochastic matrix with entries Pij(t). For a set of agents I V, with I = 
{vi, . . . ,Vk}, we define the vector of opinions xj(t) = (xy-^(t), . . . ,Xv^.(t))~^ . Given an x n matrix A 
with entries Ojj, we define the k x k matrix Aj whose entries are the a^-^-- In particular, Pi(t) is 
the matrix with entries Pviv (t)- Let us remark that Pi(t) is generally not row stochastic. However, 
if / C y is a subset of agents such that no agent in / is connected to an agent in y \ / in the graph 
G(t), then it is easy to see that 

xi(t + l) = Pi(t)xi(t) 

and Pi(t) is an aperiodic row stochastic matrix. Moreover, if I is a connected component of G(t) 
then Pi(t) is irreducible. 

The following sections are devoted to the analysis of the community structure of the network of 
agents. 

3. Asymptotic Connectivity and Agreement 

In this section, we explore the relation between communities and asymptotically connected compo- 
nents of the network. Let us remark that the set of edges E can be classified into two subsets as 
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follows: 

Ef = £ E\ 3Uj e N, Vs > Uj, (ij) ^ E{s)} 

and 

= {{iJ) e Vt G N, 3s > t, (ij) G E{s)} . 
Intuitively, an edge is in E-f if the agents i and j stop interacting with each other in finite 

time. E°^ consists of the interactions between agents that are infinitely recurrent. It is clear that 
^/ n = and E = E^ \J E°° . Also, since E and thus E^ is a finite set, there exists T G N such 
that 

(3.1) v(i,j)ei?^, Vs>r, {i,j)iE{s). 

Let us remark that the sets E^ and E'^ and the natural number T generally depend on the vector 
of initial opinions x". We define the graph = {V,E°°). 

Definition 3. Let i,j E V, we say that agents i and j are asymptotically connected if and only 
if i and j are connected in We say that they are asymptotically disconnected if they are not 

asymptotically connected. 

3.1. Asymptotic Connectivity Implies Asymptotic Agreement. 

Proposition 2. Under Assumption^ (Stochasticity) , if two agents i,j G V are asymptotically 
connected then they asymptotically agree. 

Proof. Suppose {i,j) £ E°^. From the definition of E°" there exists a strictly increasing sequence 
of non-negative integers (rfc)fcgN such that for all k £ N, {i,j) G E{Tk). Then, for all /c G N, 
\xi{Tk) — Xj{Tk)\ < Rp'^'' . Since p G (0, 1) and lim = +c« and one gets lim Xi{Tk) = lim Xj{Tk). 

On the other hand, the sequences Xi{t) and Xj{t) are convergent, which ensures that 

X* = lim Xi{t) = lim Xi{Tk) = lim Xj{Tk) = lim Xj{t) = x* 

The result in the proposition then follows from the transitivity of equality and the definition of 
asymptotic connectivity. □ 

Remark 4. The notion of asymptotic connectivity has already been considered in several works 
(including |JLM03| IBHOT051 IMorOS] ) for proving consensus in multi-agent systems. Actually, the 
previous proposition could he proved using Theorem 3 in |Mor05| . However, for the sake of self- 
containment, we preferred to provide a simpler proof of the result that uses the specificities of our 
model. 

3.2. Asymptotic Agreement Implies Asymptotic Connectivity. The converse result of Propo- 
sition[2]is much more challenging: it is clear that it cannot hold for all initial conditions. Indeed, if all 
the initial opinions x^ are identical, then it is clear that the agents asymptotically agree even though 
some of them may be asymptotically disconnected which would be the case if the graph G is not 
connected. Therefore, we shall prove that the converse result holds for almost all initial conditions. 
In this paragraph, we will need additional assumptions in order to be able to prove this result. The 
first one is the following: 

Assumption 2 (Invertibility and graph to matrix mapping). The sequence of matrices P{t) satisfy 
the following conditions: 

(a) For all t G N, P{t) is invertible. 
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(b) For all t G N, f G N, if G{t') = G{t) then P{t') = P{t). 

The first assumption is quite strong and we notice tliat it is not verified by tlie original Krause model. 
However, it can be enforced, for instance, by choosing Pii{t) > 1/2 for all i £ V, for all i G N, in 
that case P{t) is a strictly diagonally dominant matrix and therefore it is invertible. The second 
assumption states that P{t) only depends on the graph G{t), then we shall write P{t) = P{G{t)) 
where P{G') is the matrix associated to a graph G' G S{G). From the first assumption, P{G') must 
be invertible. Then, we can define for all t G N, the following set of matrices: 

(3.2) Qt = {P(Go)-^P(Gi)^^ . . . P{Gt-i)-^\ Gk G S{G), 0<k<t-l}. 

Let us remark that since S{G) is finite, the set Qt is finite: it has at most 2*^1^1/^ elements. 

We shall now prove the converse result of Proposition [2] in two different cases. 

3.2.1. Average preserving dynamics. We first assume that the opinion dynamics preserves the average 
of the opinions: 

Assumption 3 (Average preserving dynamics). For all t G N, for all j G V, ^^=iPij{t) = 1- 

This assumption simply means that the matrix P{t) is doubly stochastic. It is therefore average 
preserving: the average of x{t) is equal to the average of x{t + l). Also, if / C y is a subset of agents 
such that no agent in / is connected to an agent in y \ / in the graph G{t), it is easy to show that 
P[{G{t)) is average preserving. 

We now state the main result of the section: 

Theorem 1. // the matrices P{t) satisfy Assumptions^ (Stochasticity), [£| (Invertibility and graph 
to matrix mapping) and^ (Average preserving dynamics), for almost all vectors of initial opinions 
x^ , two agents i,j £ V asymptotically agree if and only if they are asymptotically connected. 

Proof. The if part of the theorem is a consequence of Proposition [2] To prove the only if part, let us 
define the following set 

W = {(I, J)| (I c y) A (/ ^ 0) A (J C F) A ( J ^ 0) A (/ n J = 0)} . 

Since V is a finite set, it is clear that W is finite (it has less than 2^" elements). For all (/, J) G W, 
let |/| and \J\ denote the number of elements of I and J respectively. We define the vector of M", 
cjj whose coordinates c/j^fc = 1/|/| if A; G /, cjj^k = — if A; G J, and c/j^fc = otherwise. We 
define the (n — l)-dimensional subspace of M": 

Hij= < xGM"|c/j-x = 5^Xi/|/|-^x,-/|J|=0 

Finally, let us define the subset of M": 

^° = uf U ([jQfu 

teN \(/,j)ew \QeQt 

where Qt is the set of matrices defined in (3.2). Since W is a finite set and for all t G N, Qt are finite 
sets, is a countable union of {n — l)-dimensional subspaces of M". Therefore A" has Lebesgue 
measure 0. 
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Let X G M" be a vector of initial opinions, let us assume that there exist two agents i,j G V that 
asymptotically agree but are asymptotically disconnected. Let us show that necessarily, belongs 
to the set . Let I and J denote the connected components of G°° containing i and j respectively. 
Since i and j are asymptotically disconnected, I (1 J = ^, therefore (I, J) £ W. Let T be defined 
as in equation ( |3.1[ ) (i.e. E{t) C E°°, Vt > T), since no agent in / is connected to an agent outside 
of I in G°° (and hence in G{t) for t > T), we have that for all t >T, xi{t + l) = Pi{G{t))xi{t). 
Moreover, Pj{G{t)) is average preserving. Therefore, for all t > T, the average of xi{t) is the same 
as the average of xj{T). From Proposition [2| all agents in / asymptotically agree, then the limit 
value is necessarily the average of xi{T). Therefore x* = ■ x/(T))/|/| where denote the Ill- 
dimensional vector with all entries equal to 1. A similar discussion gives that x* = (l|j| • xj{T))/\J\. 
Since i and j asymptotically agree, we have (1|/| • x/(r))/|/| = (l|j| • xj{T))/\J\. This means that 
x{T) G Hij and therefore 

xO = P(G(0))-ip(G(l))-i...P(G(T-l))-ix(r)G IJ QHij 

QeQr 

which leads to x^ G X^. □ 



Hence, in the case of average preserving dynamics, asymptotic connectivity is equivalent to asymp- 
totic agreement for almost all vectors of initial opinions. We shall now prove a similar result under 
different assumptions. 



3.2.2. Fast convergence assumption. We now replace the average preserving assumption by another 
assumption. From Proposition [T| we know that the opinion of each agent converges to its limit value 
no slower than O (/)*). This is an upper bound, numerical experiments show that in practice the 
convergence to the limit value is often slightly faster than 0(p*). This observation motivates the 
following assumption. 

Assumption 4 (Fast convergence). There exists p < p and M > such that for all i G V , for all 

t£N, 

\xi{t) -x*\< Mp\ 

Remark 5. The previous assumption always holds unless there exists i £ V such that 

limsup I log{\xi{t) - x*\) = log(/3). 

t— s>+oo I 

It should be noted that unlike Assumptions^ (Stochasticity), (Invertibility and graph to matrix 
mapping) and^ (Average preserving dynamics), it is generally not possible to check a priori whether 
Assumption\^ holds. However, numerical experiments tend to show that in practice, it does. 



The previous assumption allows us to state the following result: 

Lemma 1. Under Assumptions^ (Stochasticity) and^ (Fast convergence), there exists T' G N such 
that for all t > T' , G{t) = Moreover, = G.^. 



Proof. We shall prove the lemma by showing that there exists T' G 
E{t) ^ E°° C E'^ C E{t). Firstly, let Ti > T where T is defined as in equation ([sl]) 
t > Ti, E{t) C E 



N such that for aU t > T', 
then for all 



Secondly, let {i,j) G E°°, then agents i and j are asymptotically connected. 
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From Proposition [2| it follows that i and j asymptotically agree. Therefore, G E^^. Thirdly, let 
G E^g, then x* = x* and for ah t G N 

\Xi{t) — Xj{t)\ < \xi{t) — x*\ + \x* — x*\ + \Xj{t) — x*\ 
< \Xi{t) - X*\ + \Xj{t) - x*\ 

From Assumption |4j we have for all t G N, 

\xi{t) - Xj{t)\ < 2M/. 

Since p<p, there exists T2 G N, such that for ah t > T2, 2Mf/ < RpK Then, for all t > T2, G 
E{t). Let T' = max(ri, Ta), then for ah t > T' , E{t) = E"^ =^E<g and thus G{t) = = G^g. □ 

The previous result states that after a finite number of steps, the graph of interactions between 
agents remains always the same. Then, we can state a result similar to Theorem [TJ 

Theorem 2. Under Assumptions^ (Stochasticity), (Invertibility and graph to matrix mapping) 
and^ (Fast convergence), for almost all vectors of initial opinions x^ , two agents i,j^V asymptot- 
ically agree if and only if they are asymptotically connected. 

Proof. The if part of the theorem is a consequence of Proposition [2] To prove the only if part, let us 
define the following set associated to a spanning subgraph G' G S{G): 

W{G') = {{I, J)\ (/ C y) A (J C y) A (/ / J) A (/ G /C(G')) A (J G IC{G'))} . 

Since y is a finite set, it is clear that W(G') is finite (it has less than 2^*^ elements). Let {I, J) G 
W{G'), I = {vi, . . . ,^"1/1), J = {wi, . . . ,w\j\Y Since I and J are connected components of G', we 
have that Pi{G') and Pj{G') are aperiodic irreducible row stochastic matrices. Let ej{G') and ej{G') 
be the left Perron eigenvectors of Pi{G') and Pj{G'), respectively: 

ei{GyPi{G') = ei{Gy and ej{G') ■ l,,, = 1 

and 

ej{GYPj{G') = ej{G'y and ej(G') • = L 

We define the vector of M", c/j whose coordinates are given by cjj^^^, = ej^k if Vk G /, cij^^^ = —e j^k 
if Wk £ J and cjj^k = ii k £ V \ {I U J). We define the {n — l)-dimensional subspace of M": 

Hij{G') = {xe M"| cij{G') • X = 0} . 

Finally, let us define the subset of M": 

(3.3) ^° = U ( U ( U ( U QHjAG' 

teN \G'e5(G) \(/,J)e>V(G') \QeQt 




where Qt is the set of matrices defined in (3.2). S{G) is a finite set and for all G' G S{G), W{G') 
is a finite set. Moreover for all t G N, Qt is a finite set. Then, is a countable union of (n — 1)- 
dimensional subspaces of M"". Therefore has Lebesgue measure 0. 

Let x^ G M" be a vector of initial opinions, let us assume that there exist two agents i,j G V that 
asymptotically agree but are asymptotically disconnected. Let us show that necessarily, x^ belongs 
to the set . Let I and J denote the connected components of containing i and j respectively. 
Since i and j are asymptotically disconnected, I ^ J, therefore (/, J) G >V(G°°). Since / is a 
connected component of it follows from Lemma[l]that for all t > T', xi{t + 1) = Pi{G°°)xi{t). 
Moreover, Pj(G°°) is an aperiodic irreducible row stochastic matrix and from the Perron-Frobenius 
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Theorem (see e.g. |Sen81| ). it follows that 1 is a simple eigenvalue of and all other eigenvalues 

of Pj{G°°) have modulus strictly smaller than 1. Therefore, 

lim xKt) = (e/(G-)-x,(r'))l|7| 

and X* = ei{G^) ■ xjiT'). A similar discussion gives that x* = ej{G^) ■ xj{T'). Since i and j 
asymptotically agree, we have e/(G°°)-a;/(r') = ej{G'^)-x j{T'). This means that x{T') G Hij{G°°) 
and therefore 

x' = P{Gm-'P{G{l))-\..P{G{T'-l))-'x{T')e \J QHjj{G^) 

which leads to x^ G X^. □ 

In this section, we showed that asymptotic connectivity of agents implies asymptotic agreement and 
that under additional reasonable assumptions these are actually equivalent except for a set of vectors 
of initial opinions of Lebesgue measure 0. In other words, we can consider almost surely that the 
communities of agents correspond to the connected components of the graph G°°. Actually, we are 
confident that a similar result holds even without Assumptions |3] or |4j However, in this case, the 
set of initial opinions leading to agreement without connectivity is not necessarily a countable 
union of (n — l)-dimensional subspaces, and it can have much more complex geometrical features. 
Therefore, we leave the generalization of the results presented in this section as future work. 

In the following, under Assumptions [T](Stochasticity), [2] (Invertibility and graph to matrix mapping) 
and|4] (Fast convergence), we show that an algebraic characterization of communities can be given 
in terms of eigenvalues of the matrix associated to the graph of communities P{G^^). 

4. Algebraic Characterization of Communities 

Let G' £ <S{G), let / C F be a subset of agents such that no agent in / is connected to an agent 
in y \ / in the graph G', then Pi{G') is a row stochastic matrix. Let Xi{Pi{G')), . . . , A|/| (P/(G')) 
denote the eigenvalues of Pi{G') with Ai(P/(G')) = 1 and 

\Xi{Pi{G'))\ > \X2iPi{G'))\ >■■■> |A|c|(P/(G'))|. 

Let C € then no agent in C is connected to an agent in y \ C in the graph G<^. The following 
theorem gives a characterization of the communities in terms of the eigenvalues A2 (-Pc(G'c^)) for 

Theorem 3. Under Assumptions^ (Stochasticity), (Invertibility and graph to matrix mapping) 
and^ (Fast convergence), for almost all vectors of initial opinions x^ , for all communities C £ , 
such that |C| > 2, 

\\2{Pc{G^M < P- 

Proof. Let us consider a spanning subgraph G' € 5(G), let I = {vi, . . . with \I\ > 2, be a 

connected component of G' then Pi[G') is an aperiodic irreducible row stochastic matrix. Then, 
from the Perron- Frobenius Theorem, it follows that 1 is a simple eigenvalue of P/(G'). Therefore, 
A2(-P/(G')) / 1. Let fi{G') be a left eigenvector of Pi{G') associated to eigenvalue \2{Pi{G')). Let 
us define the vector of M", ci{G') whose coordinates are given by ci^y^{G') = fi^k{G') if Vk £ I and 
ci,k{G') = if fc S y \ /. We define the (n — l)-dimensional subspace of M": 

Hi{G') = {x e M"| ci{G') • X = 0} . 
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Finally, let us define the subset of ! 




^-'-Ul U U U 

yG'eSiG) \ieK{G'), \i\>2 \QeSt 




where Qt is the set of matrices defined in (3.2). S{G) is a finite set and for all G' G ^(G), /C(G') 
is a finite set. Moreover, for all i G N, Qt is a finite set. Then, is a countable union of {n — 1)- 
dimensional subspaces of M". Therefore has Lebesgue measure 0. 



Let be given as in equation (3.3), let G M" \ be a vector of initial opinions. Let us assume 
there is a community C G ^ with \C\ > 2, such that \X2iPc{G^^))\ > p. Let us show that necessarily, 
belongs to the set Y^. First, since x^ ^ X^, we have from the proof of Theorem [2] that C is 
a connected component of G°° = G^g. Therefore, from Lemma [!} there exists T' G N, such that 
for all t > T' , xc{t + 1) = Pc{G^g)xc{t) and Pc{G^i^) is an aperiodic irreducible row stochastic 
matrix. From the Perron-Frobenius Theorem, it follows that 1 is a simple eigenvalue of Pc'(G<^) and 
all other eigenvalues of Pc{G^g) have modulus strictly smaller than 1. Let ec(G^) be the left Perron 
eigenvector of Pc{G^g): 

eciG^^V PciG^g) = eciG^^V and ec{G^g) ■ l\c\ = 1 

Then 

lim xc{t) = Xc where x^ = {ec{G^g) ■ xc(r'))l|c|- 
Let us remark that for all t>T\ 

(4.1) xc{t + l)-xl = Pc{G<g){xc{t) - x*c). 

Let fc{G^g) be a left eigenvector of Pc{G^g) associated to eigenvalue A2(-Pc(G<g')): 

fc{G^y Pc{G^g) = X2{PciG^^))fciG^^V. 



Then, it follows from equation (4.1) that for all t>T', 

fciG^g) ■ (xcit) - x*c) = fciG^g) ■ (xc(T') - x*c)X2iPc{G^g)f-^'l 
Therefore, by the Cauchy-Schwarz inequality, we have for all t >T' 

\fc{G<g)-{xc{t)-x*c)\ 



\xc{t)-x}.\\ > 



wfciG^^n 

IfciG^) ■ {xc{T') - x}^ (,_^,) 



\\jc[G^^)\\ 

Since we assumed \X2{PciGcg))\ > p, we have for all t > T' 

(4.2) Wxcit) -xc\\> y^lG^' ' ■ 
Now, let us remark that it follows from Assumption [4] that for all t G N 

(4.3) \\xc{t)-x*c\\< y^\Mp'. 



Inequalities (4.2) and (4.3) give for all t>T' 

\fc{G^)-{xc{T')-x*c)\ , 



fciGv)\\p 



< V\C\Mp\ 
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Since p < p, the previous inequality holds for all t > T' if and only if \ fc{G'if) • {xc{T') — Xq)\ = 0. 
Therefore, fc{G^) ■ xc{T') = fclc^g) ■ {xciT') - x^) = which means that x{T') € Hc{G^g). 
Therefore, 

xO = P(G(0))-ip(G(l))-^..P(G(r-l))-ix(r')G U QHc{G^g) 

which leads to x^ G . Therefore, we have proved that for all vectors of initial opinions x^ G 
M" \ (XO U for ah communities C G ^ such that \C\ > 2, |A2(Pc(G^))| < P- We conclude by 
remarking that U is a set of Lebesgue measure 0. □ 

In this section, we showed that the community structure ^ satisfies some properties related to the 
eigenvalues of the matrix Pc{G^g), for C G ^. In the following, we use this result to address the 
problem of community detection in graphs. 



5. Application: Community Detection in Graphs 



In this section, we propose to use a model of opinion dynamics with decaying confidence to address 
the problem of community detection in graphs. 

5.1. The Community Detection Problem. In the usual sense, communities in a graph are groups 
of vertices such that the concentration of edges inside one community is high and the concentration of 
edges between communities is comparatively low. Because of the increasing need of analysis tools for 
understanding complex networks in social sciences, biology, engineering or economics, the community 
detection problem has attracted a lot of attention in the recent years. The problem of community 
detection is however not rigorously defined mathematically. One reason is that community structures 
may appear at different scales in the graph: there can be communities inside communities. Another 
reason is that communities are not necessarily disjoint and can overlap. We refer the reader to 
the excellent survey [ForlOj and the references therein for more details. Some formalizations of the 
community detection problem have been proposed in terms of optimization of quality functions such 
as modularity |NG04j or partition stability |LDB09 ). 

5.1.1. Quality functions. Modularity has been introduced in [NG04j . the modularity of a partition 
measures how well the partition reflects the community structure of a graph. More precisely, let 
G = (y, E) be an undirected graph with E symmetric and anti-reflexive. For a vertex i £ V the 
degree di of i is the number of neighbors of i in G. Let V he a partition of V. Essentially, the 
modularity Q(V) of the partition V is the proportion of edges within the classes of the partition 
minus the expected proportion of such edges, where the expected number of edges between vertex i 
and j is assumed to be didj/\E\: 



Q{V) 



1 \ ^ \ - / didj 



\E\ ^ ^ \ \E\ 



where aij are the coefficients of the adjacency matrix of G (uij = 1 if G E, aij = otherwise). 
The higher the modularity, the better the partition reflects the community structure of the graph. 
Thus, it is reasonable to formulate the community detection problem as modularity maximization. 
However, it has been shown that this optimization problem is NP-complete |BDG"'"08] . Therefore, 
approaches for community detection rely mostly on heuristic methods. In [New06j . a modularity 
optimization algorithm is proposed based on spectral relaxations. Using the eigenvectors of the 
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modularity matrix, it is possible to determine a good initial guess of the community structure of the 
graph. Then, the obtained partition is refined using local combinatorial optimization. In |BGLL08] . a 
hierarchical combinatorial approach for modularity optimization is presented. This algorithm which 
can be used for very large networks, is currently the one that obtains the partitions with highest 
modularity. 

However, modularity has the drawback that it fails to capture communities at different scales. The 
notion of partition stability |LDB09| makes it possible to overcome this limitation. Let us consider 
a continuous-time process associated with a random walk over the graph G where transitions are 
triggered by a homogeneous Poisson process. Assume that the initial distribution is the stationary 
distribution. Then, the stability at time t G M"*" of the partition V is defined as 

lev 

where p{I, t) is the probability for a walker to be in the class / initially and at time t. Stability 
measures the quality of a partition by giving a positive contribution to communities from which a 
random walker is unlikely to escape within the given time scale t. For small values of t, this gives 
more weights to small communities whereas for larger values of t, larger communities are favored. 
Thus, by searching the partitions maximizing the stability for several values of t, one can detect 
communities at several scales. 



5.1.2. Eigenvalues of the normalized Laplacian matrix. We give an alternative formulation of the 
community detection problem using a measure of connectivity of graphs given by the eigenvalues of 
their normalized Laplacian matrix. Let G = {V,E) be an undirected graph with V = {l,...,n}, 
with n > 2. For a vertex i € V, the degree di{G) of i is the number of neighbors of i in G. The 
normalized Laplacian of the graph G is the matrix L[G) given by 

' 1 if 2 = i and di{G) / 0, 



-1 



, if (i, j) e E, 

otherwise. 



Let us review some of the properties of the normalized Laplacian matrix (see e.g. |Chu97j ) . ^i{L{G)) = 
is always an eigenvalue of L{G), it is simple if and only if G is connected. All other eigenvalues 
are real and belong to the interval [0, 2]. The second smallest eigenvalue of the normalized Laplacian 
matrix is denoted /X2(L(G)). It can serve as an algebraic measure of the connectivity: fj,2{L{G)) = 
if the graph G has two distinct connected components, ^2{L{G)) = n/{n — 1) if the graph is the 
complete graph (for all i,j £V,i^ j, {i,j) G E), in the other cases fj,2{L{G)) G (0, 1]. 

Remark 6. The second smallest eigenvalue of the (non-normalized) Laplacian matrix is called al- 
gebraic connectivity of a graph. In this paper, we prefer to use the eigenvalues of the normalized 
Laplacian matrix because it is less sensitive to the size of the graph. For instance, if G is the com- 
plete graph then ii2{L{G)) = n/{n— 1) whereas its algebraic connectivity is n. 

Let P be a partition of the set of vertices V . For all / E "P, with |/| > 2, L{Gi) denotes the 
normalized Laplacian matrix of the graph G[ = {I,Ej) consisting of the set of vertices / and of 
the set of edges of G between elements of I. Let us define the following measure associated to the 
partition V 

fi2iV) = min /X2(L(G,)). 

— /eP,|/|>2 

Essentially, ^2(7^) measures the connectivity of the less connected component of G-p. 
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We now propose a new formulation of the community detection problem: 

Problem 1. Given a graph G = (V, E) and a real number 5 £ (0, 1], find a partition V of V such 
that for all I £V, such that \I\ > 2, fi2{L{Gi)) > 6 (i.e. ^(P) > 6). 

If ii2{L{G)) > 6, it is sufficient to choose the trivial partition V = {V}. If 6 > fj,2{L{G)), then 
we want to find groups of vertices that are more densely connected than the global graph. This 
coincides with the notion of community. The larger 6 the more densely connected the communities. 
This makes it possible to search for communities at different scales of the graph. Let us remark 
that Problem [T] generally has several solutions. Actually, the trivial partition V = {{1}, . . . , {n}} 
is always a solution. In the following, we show how non-trivial solutions to Problem [T] can be 
obtained using a model of opinion dynamics with decaying confidence. We evaluate the modularity 
of the partitions we obtain and compare our results to those obtained using modularity optimization 
algorithms presented in jNew061 IBGLL08] . 



5.2. Opinion Dynamics for Community Detection. Let a S (0,1/2), we consider the opinion 
dynamics with decaying confidence model given by: 



(5.1) 



Xiit + 1) 



Xi{t) + 
Xi{t) 



a 



{xjit) - Xiit)) ifNiit)^^ 



iiNi{t) =0 



where Ni{t) is given by equation (2.3). It is straightforward to check that this model is a particular 
case of the model given by equations (2.1) and (2.2) and that Assumption [T] (Stochasticity) holds. 
Moreover, since a G (0, 1/2) it follows that for all i G V, t € N, Pn{t) > 1/2. Therefore the matrix 
P{t) is strictly diagonally dominant and hence it is invertible. Also, P{t) = P{G{t)), where for a 
subgraph G', P{G') = Id — aQ{G') where Id is the identity matrix and 



(5.2) 



Q^,iG') 



d,{G') 




if i = j and di{G') / 0, 
if {iJ)GE', 
otherwise. 



where di{G') denotes the degree of i in the graph G' . Therefore, Assumption [2] (Invertibility and 
graph to matrix mapping) holds as well. Let us remark that the matrix P{t) is generally not average 
preserving and therefore Assumption [3] does not hold. 

Before stating the main result of this section, we need to prove the following lemma : 

Lemma 2. Let V be a partition ofV, I £V such that \I\ > 2. Then, A is an eigenvalue of PiiG-p) 
if and only if = (1 — A)/a is an eigenvalue of L{Gi). 



Proof. First, let us remark that Pi{G-p) = Id — aQ{Gj) where Q{Gi) is defined as in equation (5.2) 
Then, let us introduce the matrices R{Gi) and D{Gi) defined by 

^ if i = i and di(G7) / 0, 



Rij{Gi 



if (i,j) G Ei, 



, otherwise. 



and 



DijiGi 



VMGi) if i = j, 
otherwise. 
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Let US remark that L(G/) = D{Gi)R{Gi) and Q{Gi) = R{Gi)D{Gi). It follows that L{Gi) and 
Q{Gi) have the same eigenvalues. The stated result is then obtained from the fact that the matrix 
Q{Gi) = {Id-Pi{Gv))/a. □ 

We now state the main result of the section which is a direct consequence of Theorem[3]and Lemma[2j 
Corollary 1. Let p = 1 — a6, under Assumption^ (Fast convergence), for almost all vecto rs of 



initial opinions , the set of communities 'tf obtained by the opinion dynamics model (5.1) is a 
solution to ProblemUl 

5.3. Examples. In this section, we propose to evaluate experimentally the validity of our approach 
on three benchmarks taken from [New06] . 

5.3.1. Zachary karate club. We propose to evaluate our approach on a standard benchmark for 
community detection: the karate club network initially studied by Zachary in |Zac73j . This is a 
social network with 34 agents shown on the top left part of Figure [T] The original study shows the 
existence of two communities represented on the figure by squares and triangles. 



We propose to use our opinion dynamics model ( |5.1[ ) to uncover the community structure of this 
network. We chose 4 different values for S and 2 different values for parameters R and a. The 
parameter p was chosen according to Corollary [TJ p = 1 — a6. For each combination of parameter 
value, the model was simulated for 1000 different vectors of initial opinions chosen randomly in 
[0, 1]^^. Simulations were performed as long as enabled by floating point arithmetics. 

The experimental results are reported in Table [TJ For each combination of parameter value, we 
indicate the partitions in communities that are the most frequently obtained after running the opinion 
dynamics model. For each partition we give the number of communities in the partition, the 
measure ^2(^)5 this value being greater than 6 indicates that Problem [l] has been solved. We 
computed the modularity Q(^) in order to evaluate the quality of the obtained partition. We also 
indicate the number of times that each partition occurred over the 1000 simulations of the opinion 
dynamics model. 

We can check in Table [T] that all the partitions are solutions of Problem [T} Let us remark that in 
general the computed partition depends on the initial vector of opinions, this is the case for 5 = 0.3 
and 5 = 0.4. Also, changing the parameters R and a seems to have some effect on the probability of 
obtaining a given partition. For instance, for 5 = 0.3, the probabilities of obtaining one partition are 
significantly different for R = 1 and R = 10. Also, for 6 = 0.4, the probabilities are slightly different 
for a = 0.1 and a = 0.2. 

However, it is interesting to note that the partitions that are obtained for the same value of parameter 
6 have modularities of the same order of magnitude which seems to show that these are of compa- 
rable quality. The partition with maximal modularity is obtained for 5 = 0.4, it is a partition in 4 
communities with modularity 0.417. As a comparison, algorithms |New06l IBGLL08] obtain a parti- 
tion in 4 communities with modularity 0.419. This shows that our approach not only allows to solve 
Problem [1] but also furnishes partitions with a good modularity which might seem surprising given 
the fact that our approach, contrarily to |New061 IBGLLOSj does not try to maximize modularity. 

In Figure [TJ we represented the graphs of communities G-^ that are the most frequently obtained for 
R = 1, a = 0.1 and the different values of 6. It is interesting to remark that for 6 = 0.2 we almost 
obtained the communities that were reported in the original study |Zac73j . Only one agent has been 
classified differently. One may argue that this agent has originally 4 neighbors in each community 
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6 


1^1 






Occurences 
K = 1, a = 0.1 


Occurences 
i? = 10, a = 0.1 


Occurences 
K = 1, a = 0.2 


Occurences 
i? = 10, a = 0.2 


0.1 


1 


0.132 





1000 


1000 


1000 


1000 


0.2 


2 


0.250 


0.360 


1000 


999 


1000 


999 


0.3 

0.3 


CO CO 


0.334 
0.363 


0.399 
0.374 


691 
283 


105 
891 


679 
298 


63 
937 


0.4 
0.4 


4 
5 


0.566 
0.566 


0.417 
0.402 


924 
15 


994 
6 


884 
54 


897 
98 



Table 1. Properties of the partitions of the karate club network obtained by the 
opinion dynamics model (1000 different vectors of initial opinions for each combination 
of parameter values). 




Figure 1 . Graphs G^g- for the most frequently obtained partition of the karate club 
network for = 1, a = 0.1 and 5 = 0.1 (top left), 5 = 0.2 (top right), 6 = 0.3 (bottom 
left), 6 = 0.4 (bottom right). 



so it could be classified in one or the other. It is also interesting to see that our approach allows 
us to search for communities at different scales of the graph. When 5 increases, the communities 
become smaller but more densely connected. This is corroborated by computing the stability of these 
partitions (see Figure [2]) . We can see that the partition with maximal stability changes according to 
time-scale t: for small values of t the partition in 4 communities is better, for intermediate values 
of t the partition in 3 communities has the largest stability, for large values of t the partition in 2 
communities maximizes the stability. 
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0.7r 



8=0.2 

S=0.3 

5=0.4 




10 10 10 

Figure 2. Stability of the partitions presented in Figure [T] 



5.3.2. Books on American Politics. We propose to use our approach on an example consisting 
of a network of 105 books on politics |New06| . initially compiled by V. Krebs (unpublished, see 
www. orgnet . com). In this network, each vertex represents a book on American politics bought from 
Amazon.com. An edge between two vertices means that these books are frequently purchased by 
the same buyer. The network is presented on the top left part of Figure |3] where the shape of the 
vertices represent the political alignement of the book (liberal, conservative, centrist). 



We used our opinion dynamics model (5.1 ) to uncover the community structure of this network. We 
chose 3 different values for 5 and 2 different values for parameters R and a. The parameter p was 
chosen according to Corollary [T] p = l — a6. For each combination of parameter value, the model was 
simulated for 1000 different vectors of initial opinions chosen randomly in [0, 1]^^^. Simulations were 
performed as long as enabled by floating point arithmetics. The experimental results are reported 
in Table d 



6 


1^1 




Q(^) 


Occurences 

/J = 1, Q- = 0.1 


Occurences 
_R = 10, a = 0.1 


Occurences 
R = 1, o = 0.2 


Occurences 
_R = 10, « = 0.2 


0.1 
0.1 


2 
2 


0.134 
0.129 


0.457 
0.457 


980 

20 


1000 



640 
360 


581 
419 


0.15 
0.15 


CO CO 


0.182 
0.187 


0.499 
0.494 


898 
102 


1000 



905 
95 


1000 



0.2 
0.2 
0.2 


4 
4 
4 


0.269 
0.266 
0.269 


0.523 
0.512 
0.520 


678 
218 
49 


1000 




673 
207 
72 


1000 





Table 2. Properties of the partitions of the books network obtained by the opinion 
dynamics model (1000 different vectors of initial opinions for each combination of 
parameter values). 
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Figure 3. Graphs for the most frequently obtained partition of the books net- 
work for = 1, a = 0.1: initial graph (top left), 5 = 0.1 (top right), 5 = 0.15 
(bottom left), 5 = 0.2 (bottom right). Shapes represent political alignment of the 
books: circles are liberal, squares are conservative, triangles are centrist. 

Let us remark that the computed partitions are solutions to the Problem [T| Also, for the same 
value of parameter 5, the modularity is very similar for all partitions. Actually, all the partitions 
obtained for the same value of 6 are almost the same. As in the previous example, we can see that 
the choice of parameters R and a affects the probability of obtaining a given partition. The partition 
with maximal modularity is obtained for 5 = 0.2, it is a partition in 4 communities with modularity 
0.523. As a comparison, algorithms |New06j and |BGLL08] obtain partitions in 4 communities with 
modularity 0.526 and 0.527, respectively. As we can see, our partition has a modularity that is quite 
close from those obtained by these algorithms. 

In Figure |3j we represented the graphs of communities that are the most frequently obtained for 
the different values of 6. Let us remark that even though the information on the political alignment 
of the books is not used by the algorithm, our approach allows to uncover this information. Indeed, 
for 6 = 0.1, we obtain 2 communities that are essentially liberal and conservative. For 6 = 0.2, we 
then obtain 4 communities: liberal, conservative, centrist-liberal, centrist-conservative. 

In Figure |4j we represented the stability of the partitions shown in Figure m As in the previous 
example, we can see that the partition with maximal stability changes according to time-scale t 
which shows that our approach makes it possible to detect community at several scales using different 
values of parameter 5. 



OPINION DYNAMICS WITH DECAYING CONFIDENCE 



19 



0.7 



0.6 



8=0.1 

8=0.15 

8=0.2 



O 
en 




OL 

10" 



-1 



Figure 4. Stability of the partitions presented in Figure [3| 



5.3.3. Political blogs. The last example we consider consists of a significantly larger network of 1222 
political blogs |AG05j . In this network, an edge between two vertices means that one of the corre- 
sponding blogs contained a hyperlink to the other on its front page. We also have the information 
about the political alignment of each blog based on content: 636 are conservative, 586 are liberal. 

The two previous examples show that the modularity of the obtained partitions does not depend 
much on the parameters R and a or on the vector of initial opinions. For this reason, we decided 
to apply our opinion dynamics model with parameters R = 1 and a = 0.1. We used 17 values of 
6 between 0.05 and 0.75. The parameter p was chosen according to Corollary [T] p = 1 — a6. For 
each value of 6, the model was simulated only once for a vector of initial opinion chosen randomly 
m 

[0,1] 1222^ Simulations were performed as long as enabled by floating point arithmetics. 

The partition with maximal modularity was obtained for 6 = 0.4. It is a partition in 12 communities 
with modularity 0.426. There are 2 main communities: one with 653 blogs, from which 94% are con- 
servative, and one with 541 blogs, from which 98% are liberal. The 28 remaining blogs are distributed 
in 10 tiny communities. When we progressively increase 6, we can see that the size of the two large 
communities reduces moderately but progressively until 6 = 0.65 where the conservative community 
splits into several smaller communities, the largest one containing 40 blogs. The liberal community 
remains until 5 = 0.725 where it splits into smaller communities, the largest one containing 54 blogs. 

As a comparison, algorithm |New06| obtains a partition in 2 communities with modularity 0.426 
whereas algorithm |BGLL08] obtains a partition in 9 communities with modularity 0.427. As we can 
see, the partition we obtain is very acceptable in terms of modularity. 

In Table |3} we give a comparative summary of the modularity of the partition obtained for the three 
examples by our approach and by the algorithms presented in [New061 IBGLLOSj . Though slightly 
smaller, the modularity of the partition we obtain is comparable to that of other partitions which 
is actually surprising since our approach, contrarily to jNew061 IBGLL08] does not try to maximize 
modularity. 
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Network 


Karate 


Books 


Blogs 


Number of nodes 


34 


105 


1222 


This article 


0.417 


0.523 


0.426 


|New06j 


0.419 


0.526 


0.426 


IBGLL08I 


0.419 


0.527 


0.427 



Table 3. Modularity of the partitions obtained by the approach presented in this 
paper and by the algorithms presented in |New061 IBGLLOSj for the three examples 
considered in this paper. 



6. Conclusion and Future Work 

In this paper, we introduced and analyzed a model of opinion dynamics with decaying confidence 
where agents may only reach local agreements organizing themselves in communities. Under suit- 
able assumptions, we have shown that these communities correspond to asymptotically connected 
components of the network. We have also provided an algebraic characterization of communities in 
terms of eigenvalues of the matrix defining the collective dynamics. To complete the analysis of our 
model, future work should focus on relaxing Assumption [4] by studying the model behavior when 
there is an agent i that approaches its limit value at a rate exactly p: 

limsup \ \og{\xi{t) - x*\) = log(p). 

In the last part of the paper, we have applied our opinion dynamics model to address the problem of 
community detection in graphs. We believe that this new approach offers an appealing interpretation 
of community detection: communities are sets of agents that succeed to reach an agreement under 
some convergence rate constraint. We have shown on three examples that this approach is not only 
appealing but is also effective. In the future, we shall work on a distributed implementation of 
our approach. Let us remark that this should be feasible since our approach is by nature based 
on distributed computations. Then, we shall use our approach to analyze a number of networks 
including large scale networks. 
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